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Abstract 

We introduce an arrival time operator which is self-adjoint and, unlike previously proposed arrival 
time operators, has a close link to simple measurement models. Its spectrum leads to an arrival 
time distribution which is a variant of the Kijowski distribution (a re-ordering of the current) in the 
large momentum regime but is proportional to the kinetic energy density in the small momentum 
regime, in agreement with measurement models. A brief derivation of the latter distribution is 
given. We make some simple observations about the physical reasons for self-adjointness, or its 
absence, in both arrival time operators and the momentum operator on the half-line and we also 
compare our operator with the dwell time operator. 
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I. INTRODUCTION 


The arrival time problem in quantum mechanics is the question of determining the prob¬ 
ability that an incoming wave packet, for a free particle, arrives at the origin in a given time 
interval [1-5]. Classically, for a particle with initial position x and momentum p, the arrival 


time is given by the quantity 


r = 


mx 


P 


( 1 ) 


The quantum problem is most simply solved using the spectrum of an operator corresponding 
to this quantity, such as that first studied by Aharonov and Bohrn [6], 
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A heuristic result due to Pauli |7j (significantly updated by Galapon {§,9]) indicates that 
an object such as this, which is conjugate to a Hamiltonian with a semi-bounded spec¬ 
trum, cannot be self-adjoint. Indeed we find that its eigenstates, which in the momentum 
representation (with x —> ihd/dp ) are given by 


0r(p) = 


\p\ v 
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jJP—n 

g 2 mfi 


( 3 ) 

from 


are complete but not orthogonal. There is a POVM associated with these states jioj] 
which an arrival time distribution can be constructed and it coincides with that postulated 
by Kijowski 11], namely 


n k ( t ) = | (V#t)P 


= — (^t\WKx)\p\ 2 Ht) 
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( 4 ) 


(where n(r)dr is the probability of arriving at the origin between r and r + dr ), for which 


there is some experimental evidence 


12]. This is related by a simple operator re-ordering 


to the quantum-mechanical current at the origin, (J(f)), the result expected on classical 
grounds, where the current operator is given by 


J(t) = 2^ (*)) + S{x{t))p ), 


( 5 ) 


with x(t) = x +pt/m. This picture, the standard one, is nicely summarized in Ref. 13] and 


some aeve 


Refs. 


opments of it and explorations of the underlying mathematics are described in 
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The purpose of the present paper is to make two contributions to the standard pic¬ 
ture presented above. The first is to discuss three simple self-adjoint modifications of the 
Aharonov-Bohm operator, discuss the relationship with the momentum operator on the half¬ 
line, and identify the underlying physically intuitive reasons why some of these operators are 
self-adjoint and some not. The second and main result is to present a new self-adjoint arrival 
time operator which has a much closer link to models of measurement than any previously 
studied operators. 


II. SELF-ADJOINT ARRIVAL TIME OPERATORS. 


The lack of self-adjointness of Eq. (l2jl is not necessarily a problem but nevertheless a 
number of efforts have been made to restore it. Here, we take a simple approach and note 
that self-adjointness may be achieved by a number of simple modifications of the states 
Eq. (J3]). The states 

/ -M_V e i<P)^K\ (6) 
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where e(p) is the sign function, are orthogonal and complete and so are eigenstates of a self- 
adjoint operator. This operator, first considered by Kijowski [ll] and subsequently explored 


at length by Delgado and Muga 


16] . may be written 
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( 7 ) 


and is a quantization of the classical expression —mx/\p\. It is also usefully written in the 
representation using the pseudo-energy, £ = p\p\/2m, where we have 

d 


T, 


KDM 


= —ih 


di 


( 8 ) 


(which acts on states $(£) = (m/p)?4>(p)), and is self-adjoint since £ takes an infinite range. 
This is in contrast to the Aharonov-Bohm operator which, in the energy representation, has 
the form 

Tab = ( ® ( -ih-^ ) (9) 


dE J \ dE / 

where the two parts of the direct sum refer to the positive and negative momentum sectors. 
Its lack of self-adjointness is due to the fact that E > 0, as is frequently noted. (See for 
example, Ref. jjj]). 
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A second modification of the Aharonov-Bohm operator is to superpose opposite val¬ 
ues of r in Eq.(J3J) and then note that the subsequent states, which are proportional to 
\p\% sm(p 2 r/2mh) are orthogonal and are the eigenstates of the self-adjoint operator 


-MI 


= \TI b , 


( 10 ) 


considered by de la Madrid and Isidro 18[. This is a quantization of m|x|/|p|. A third mod¬ 
ification is to note that the orthogonality of the states \p\% sin(p 2 T /2mh) is not affected by 
restriction to positive or negative momenta so we may consider these two sectors separately 
and as a consequence the operator 


T 3 = 0(p)T mi 0(p) - 0(-p)T MI 9(-p), 


( 11 ) 


is self-adjoint. This operator, which does not seem to have been noted previously, is a 
quantization of the classical expression m\x\/p. 

These three examples all side-step the Pauli theorem since they do not have canonical 
commutation with the Hamiltonian. Furthermore, they all give probability distributions 
which are simple variants on the Kijowski distribution, the expected result, as is easily 
deduced from their eigenstates. 

From these three examples we make the following simple observation. The Aharonov- 
Bohm operator arises from the quantization of the classical expression — mx/p and is not 
self-adjoint. However, quantizing any of the three classical expressions —mx/\p\, m|x|/|p| or 
m\x\/p leads to a self-adjoint operator. Hence, self-adjoint modifications of the Aharonov- 
Bohm operator are easily obtained by relinquishing just one or two bits of information, 
namely the sign of x, or p, or the signs of both. The relinquished information is essentially 
the specification of whether the particle is incoming (x and p with opposite signs) or outgoing 
(x and p with the same sign). Of course in practice we are usually interested in the arrival 
time for a given state, for which this information is already known, at least semiclassically, 
so from this point of view the difference between the Aharonov-Bohm operator and its 
self-adjoint variants may not be important. Nevertheless it is of interest to uncover the 
underlying origins of self-adjointness or its absence and the above properties give some 
useful clues. 

Physically speaking, self-adjointness or its absence are about precision. A self-adjoint 
operator has orthonormal eigenstates and an associated projection operator onto a range 


4 




of its spectrum. Projections onto different ranges have zero overlap. An operator that is 
not self-adjoint has non-orthogonal eigenstates and has at best a POVM onto a range of 
its spectrum. Two POVMs localizing onto different ranges will have a small overlap which 
means there is intrinsic imprecision in the specification the ranges they localize onto. To 
understand the origin of the lack of self-adjointness in the Aharonov-Bohm operator we 
would like to find a physically intuitive understanding of where this imprecision comes from. 


III. THE MOMENTUM OPERATOR ON THE HALF-LINE 

To understand the above issue, we turn to the frequently-studied situation of the mo¬ 
mentum operator on the half-line x > 0 


19 


-21J. There, the momentum operator cannot be 
made self-adjoint since it generates translations into negative x. However, p 2 can be made 
self-adjoint, with suitable boundary conditions, and therefore, by the spectral theorem, \p\ 
can be made self-adjoint. Hence just by relinquishing information about the sign of p a self- 
adjoint operator is obtained. The obstruction to self-adjointness on the half-line therefore 
lies in the sign function of p. Differently put, the problem is that the operator 8(x)9(p)6(x) 
cannot be made self-adjoint. For similar reasons, we also note that the position operator on 
the positive momentum sector cannot be self-adjoint. A POVM for the momentum operator 


on the half-line may be constructed, but this is not directly relevant to what we do here 


21 ]- 


We propose that there is a simple physical way of understanding the underlying impre¬ 
cision linked to this lack of self-adjointness. Suppose we tried to measure the momentum. 
Let us therefore consider a simple measurement model of momentum on the half-line x > 0 
using sequential position measurements, from which information about momentum can be 
deduced. Similar approaches to calculating the time-of-flight momentum have been given 
elsewhere 22] and we make use of these results, but adapted to the case of propagation in 


the region x > 0. We suppose we have an initial incoming state i/j at time t 0 consisting of 
spatially very broad gaussian, close to a plane wave, of momentum p 0 < 0 and we ask if it 
passes through a spatial region [x\ — A,fh + A] in x > 0 at time t\ and at a later time t 2 
through a spatial region [x 2 — A, x 2 + A]. The probability for these two measured results is 


p{x i, h, x 2 , t 2 ) = (01 g\h, t 0 )P Sl g*(t 2 , ti)Px 2 g(t 2 , tO-Px^ti, t 0 )\ip) 


( 12 ) 
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where 


ex+A 


Prr = 


dx\x)(x\ 


(13) 


' x—A 


is a projector onto the range [x — A, x + A] and g(ti, t 0 ) denotes the propagator in the region 
x > 0. The precise form of the propagator depends on the boundary conditions on the states 
imposed at x — 0. There is a one-parameter family that leads to a self-adjoint Hamiltonian, 
of the form 

-0(0) + iatp'ifi) — 0 (14) 

where a is real [h]]. For simplicity, we focus on the simplest case a = 0 in which the states 
are required to vanish at x = 0 and in this case the propagator in x > 0 is given by the 
simple method of images expression 


g(xi, ti\xo, t 0 ) = 9(x 1 )d(x 0 ) 
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exp 


2nh(ti — t 0 ) 
im(x i — xo) 2 


exp 


im(x i + xqY 
2h(t l - 1 0 ) 


(' 15 ') 

2/z(ti — t 0 ) ) * V 2h(ti — to) )\ 

Semiclassically, the first term clearly corresponds to the direct classical path between the 
initial and final spacetime points and the second term corresponds to a reflected path (or 
equivalently, to the path coming from the image point <—xo). Other choices of boundary 
conditions yield different propagators but all have the same feature of including both direct 
and reflected paths. 

To determine the momentum from the probability Eq. ([T2]) we consider the conditional 
probability of finding the particle near x 2 at time t 2 , given that it was near x\ at time ti, 
which is given by 


P{x2,t 2 \xi,tl) = 


p(x 1 ,t 1 ,X 2 ,t 2 ) 
p(xi,h) 


(16) 


By comparing with earlier cases for the whole real line 


22], it is reasonably easy to see 


what happens without doing any further calculation. This probability distribution will have 
two peaks, one at the expected value for an incoming state with momentum p 0 , namely 
x 2 = xi + (po/m)(t 2 — ti). The other peak is at x 2 = x\ — ( po/m){t 2 — ti) corresponding 
to the parts of the incoming state that already reached the origin prior to t\ and been 
reflected back. Hence the measured momentum probability will be concentrated about both 
p 0 and — po without a possibility of distinguishing between the two. The magnitude of 
the momentum is measured unambiguously but the sign of the momentum is intrinsically 
ambiguous, due to reflection. 
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Of course the case described above is the extreme case of an incoming plane wave with 
no spatial localization so the state arrives at the first measurement at a very wide spread of 
times. The opposite extreme is a wave packet tightly peaked in both position and momentum 
which arrives at the first measurement at a precise time and the probability Eq. (1121) is 
strongly peaked about the direct momentum p 0 with negligble peak about the reflected 
one. But for a more general state, there always will be some peaking about the reflected 
momentum and hence ambiguity in the momentum measurement. This is the physically 
intuitive reason underlying the imprecision in the sign of the momentum operator on the 
half-line associated with its lack of self-adjointness. 

This intuitive link between reflection and lack of self-adjointness is further substantiated 
in another closely related but different case, namely the momentum operator in a finite 
spatial interval [a, b\. In this case the momentum operator is in fact self-adjoint, for suitable 
periodic boundary conditions on the state, the simplest of which is 


</>(“) = V’(U 


( 17 ) 


(see for example Ref. |l9|). One can again ask how this reconciles with measurement of mo¬ 
mentum in this interval, using sequential position measurements - is there any imprecision 
due to reflection effects? The point here is that the boundary conditions mean the config¬ 
uration space has the topology of the cylinder and the semiclassical interpretation of the 
propagator is that it consists of straight line paths on the cylinder. There is no reflection at 
the boundaries, as there is in the half-line case. This is consistent with self-adjointness. 

Given these features of the momentum on the half-line, we may now draw a direct con¬ 
nection with the Aharonov-Bohm operator. In exploring the quantization of the classical 
arrival time —mx/p, it is not at all obvious at the classical level why this object leads to 
an operator which is not self-adjoint. To explore this, one might take as a simple starting 
point the simple expression 9{t + mx/p), which may be used to address the simple question 
of whether the arrival time is greater or less than a given time t. The arrival time is then 
constructed using 



(18) 


Quantizing 9(t + mx/p) directly is problematic since —mx/p does not turn into a self-adjoint 
operator so it is natural to consider alternative but equivalent classical expressions obtained 
by scaling out p to relate it to the expression 9(x+pt/m), which clearly can be turned into a 
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self-adjoint operator. (The latter expression is also a natural starting point for measurement 
inspired models as we shall see.) However, the sign of p is important in this scaling and we 
in fact obtain 

9 (* + = 9 (p) 9 ( x + + 9 (-^ e (“* “ ( 19 ) 

This is just 9(x)9(p) + 9(—x)9(— p), shifted along the classical equations of motion, precisely 
the object that cannot be made self-adjoint as we have seen. For example, with a particular 
choice of operator ordering the quantization of the right-hand side is 


exp 



[9(x)9(p)9(x) — 9{—x)9(— p)9 (—£)] exp 



( 20 ) 


A different operator ordering would lead to the expressions of the form 9(p)9(x)9(p) which 
are also not self-adjoint, as noted already. Hence, the quantization of the classical arrival 
time expression —mx/p boils down to the expression 9(x)9(p)9(x), precisely the underlying 
non-self-adjoint object encountered in the momentum operator on the half line. The same 
remarks concerning reflection therefore apply. 

By contrast, for the classical expression — mx/\p\, similar steps yield 


9 t + 


mx 

\P\ 


= 9 \ x + 


\p\t 


m 


( 21 ) 


The right-hand side, when quantized, is simply 9(x) unitarily evolved in time with the 
pseudo-Hamiltonian £ = p\p\/2m, so is self-adjoint. 

In summary, the key difference between the Aharonov-Bohm operator and the three self- 
adjoint variants of it is that the Aharonov-Bohm operator involves specification of whether 
the states are incoming or outgoing. This involves specifying the signs of the momenta on 
the two half-lines x > 0 and x < 0. Restriction to the half-line in quantum mechanics 
creates quantum-mechanical reflection which renders the sign of the momentum imprecise. 
This is the physical reason why we would expect that the Aharonov-Bohm operator is not 
self-adjoint. 


IV. DIFFICULTIES WITH THE STANDARD ARRIVAL TIME OPERATORS 

Aside from the self-adjointness issue, all four of the operators discussed in Section II 
suffer from a number of problems. Firstly, their relation to actual measurements, or at least 
to simple measurement models is not obvious since there is no physical system that couples 



to the Aharonov-Bohm operator or any of the above variants of it. Secondly all the above 
arrival time operators yield an arrival time distribution which is the Kijowski distribution or 
variants thereof for all ranges of momenta. Simple measurement models (such as those based 


on a complex potential [12 


.23] or stopwatch 


24 


25J) agree with this for large momenta. 


But crucially, for small momenta, and more specifically for arrival time measurements more 
precise than hm/p 2 (the energy time 26,(27]), reflection off the detector becomes significant 


and measurement models typically yield a distribution proportional to the kinetic energy 
density 

n(r) = N(il> T \p6{x)p\il> T ), (22) 

where N is a model-dependent normalization factor 0 , 0 - 


-25 


28]. This behaviour does 


not arise in any of the standard arrival time operators at small momentum. 

Thirdly, even after normalizable states are constructed from the eigenstates of Tab and 
its three modifications (by superposing over narrow ranges of r, as described in Ref. 25]), 
they all go like p 3 for small p, which means that they have infinite (Arc) 2 , so have poor 
spatial localization properties, contrary to the intuitive notion of what an arrival time state 
should look like. The last two problems are issues around small momentum behaviour and 
could be solved by an arrival time operator whose eigenstates go like p for small p. The 
main result of the remainder of this paper, and the solution to these three problems, is the 
construction of an arrival time operator directly inspired by simple measurement models. 

We briefly mention the work of Grot, Rovelli and Tate, who removed the singularity 
at p = 0 in Tab by a somewhat artificial regularization procedure 29]. This produced a 
self-adjoint operator but the low momentum behaviour produces a spatial spread even more 


severe than the examples above 


25], reiterating the need for a physically motivated handling 


of the low momentum regime. See also Ref. 
operators. 


30] for other criticisms of traditional arrival time 


V. MEASUREMENT MODELS 

To motivate our proposed new arrival time operator and also to substantiate the arrival 
time distribution in the low momentum regime Eq. (j22l) we consider two measurement models 
for arrival time. 

A simple measurement model for the probability of a particle crossing the origin in a 
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given time interval [0,r] involves spatial measurements onto x > 0 and x < 0 described by 
projectors P = 9(x) and P = 9(—x), and we simply check to see if the particle is on opposite 
sides of the origin at the initial and final times. The probability for crossing is then 


p(0, r) = (iJ>\PP(t)P\iJ>) + (^\PP(r)P\^). (23) 

Since dP(t)/dt = J{t), where J(t) is the current operator Eq. (15|) . this may be rewritten 


p ( 0 , t )= / dt(4>\pj(t)p\ip)- / dt(ip\pj(t)p\il>), 


(24) 


thus indicating the appearance of the current operator in arrival time probabilities derived 
from a measurement model (as it does in models, e.g. Ref. [2d|). This is not in fact the 
Kijowski distribution (integrated over a time range) Eq.flU). However, estimates involving 
gaussian states show that it is very close either for states strongly peaked in momentum, or 
for time intervals large compared to the energy time, hm/p 2 . 

We may use Eq. ([2Tj) to see the origin of the low momentum regime formula Eq. (l22l) as 
follows. The formula Eq. f|22j) is the appropriate one in the case of very frequent measurement 
in which there is significant reflection off the detector. To model this, we suppose that prior 
to t — 0, the system has been subjected to frequent projections onto x < 0, so that the state 
at t — 0 has the form 

\^)=P(t n )---P(t 2 )P(t 1 )\ ( p 0 ) (25) 

for a sequence of times 0 > t n > ■ ■ ■ > £2 > H and for some initial state |0o) hr the 
distant past. If the interval between these projections is sufficiently small (in comparsion to 
the energy time hm/p 2 ) the evolution in Eq. (l25j) will in fact be well-approximated by the 
restricted propagator for the half-line x < 0, given by Eq. (lT5j) adapted to the case x < 0, 


so there will be reflection off the origin 


26]. At t — 0 the state will therefore have the 


approximate form ^(x) = 9(—x)(4>{x) — (p(—x)), for some 4>(x), so ijj(x) is zero at x = 0 but 
has non-zero derivative. The arrival time distribution is then 


H(t) = (^|J(r)|^) 

= --[r(0,r)n0,r)-mrm0,r)}, (26) 

where the dash denotes spatial derivative. This is zero at r = 0 so we have to expand the 
wave function for small r. We make use of the free particle propagator and write 

^(0 ,t)= j dy ( 7 ^^:) 2 (27) 
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where we have used the fact that ip(x > 0, 0) = 0. Making the change of variables y = zr 2 , 
expanding 0 (zt 2 ,O), and recalling that "0(0,0) = 0, we thus obtain 

0(°,^) ~ (^^) 2 0 , (°,°)/ dzze tt. (28) 

Evaluating the integral we then find 

W)“0r(-)b 4 1^(0,0)| 2 . (29) 

2 tt 2 \mj 

This is of the desired form, Eq. 022D since 


(0|p<5(£)0|0) = h 2 |0'(O,i 


(30) 


(A similar derivation was given in Ref.|31| with different aims). Hence the low momentum 
regime result arises simply because the very frequent measurement causes the wave function 
to vanish at the origin so it is necessary to expand the average current around this for small 
times. 

A second simple model for measuring the arrival time, considered by numerous authors 


24 


25], consists of a stopwatch - a system with coordinate y and zero Hamiltonian which 


couples to the particle through the interaction p y 9(—x). Classically, for a particle approach¬ 
ing from the left it therefore causes a shift in the stopwatch variable y for the entire time 
the particle is in x < 0, stopping when the particle reaches the origin, with final value 


y(T)-y( 0 ) = 


dt 9 ( —x — —t 
m 


(31) 


where T is taken to be very large. This is easily seen to be equal classically to —mx/p 
for p > 0. In the quantum-mechanical analysis of this system we can see the form of the 
interaction between particle and stopwatch via the S'-matrix, 

r T 


S = Texp ( — — A I dt p y (t)9 ( —x — —t 


(32) 


where T is the usual time ordering operator and A is the coupling constant. The stopwatch 
has zero free Hamiltonian so p y (t) is independent of t. For weak couplings we therefore 
see that the combination of the particle variables that the stopwatch “sees” is precisely 
the quantization of the right-hand side of Eq. (1311) . It does not couple to anything like the 
Aharonov-Bohm operator. 


11 








It is also of interest to rewrite the right-hand side of Eq. (I3TT) to indicate its connection 
to the current. We have, 


[ dt 9 (—x — — dt tJ(t ) (33) 

Jo v m 

where J(t) = ( p/m)5(x + pt/m ) is the classical current and we have dropped a boundary 


L 



(33) 


term since 9(—x — pt/m) is zero for large t with positive p. We could also contemplate a 
similar integration by parts for the corresponding quantum operators but there it is less 
obvious that the boundary term may be dropped. Nevertheless, both sides of Eq. (133|) give 
two alternative starting points for the construction of arrival time operators, classically 
equivalent to —mx/p, but potentially different in the quantum case and arguably more 
relevant to measurement models. The expression on the left-hand side is difficult to handle 
as a quantum operator so we focus on a version of the right-hand side in what follows. 

VI. NEW OPERATOR AND ITS PROPERTIES 

Motivated by the above observations and in particular by the key role the current plays 
in measurement models, we look for a new time operator defined in terms of the current 
operator J(t), Eq. ()5il . We begin by noting the classical result indicated already, namely 



(34) 


where J{t) is the classical current, here valid for all values of x and p. Quantization of 
the left-hand side yields the KDM operator Eq. (I7)) . Here, however, we instead take the 
right-hand side as the starting point for quantization and this leads to a different operator, 
namely 



(35) 


whose properties we study. To be clear, there are two motivations for examining this expres¬ 
sion. Firstly, it has a potentially closer connection to measurements than the KDM operator, 
as indicated. But secondly, and independently of this, it may be simply regarded as an ex¬ 
ploration of the consequences in the quantum theory of taking a different but classically 
equivalent starting point. 

To evaluate the integral in Eq. fl35l) we first note the operator identity 



( 36 ) 
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The remaining time integral may be evaluated by sandwiching between momentum states. 
We have 


dt (pi|<S ( x+ — ) | p 2 ) = {pi\S(x)\p 2 ) / dt exp 
rnj J -oo 

2 J2\ 




= 2m S(p 1 -p 2 ) 


m 


= |—r (£(pi -P2) + 5 (pi +P2)) 
\Pi\ 

(1 + R) 

— \P l| pq- \P2/ 


(37) 


where we have introduced the reflection operator R , defined by R\p) — \ — p). We thus find 
that our time operator may now be written 


T ~ + R> + ¥\ [1 + R)i 
It may be written in terms of the KDM operator Eq.Q as 

- - ihrri - 

T = T kdm + R 

2p\p\ 


(38) 


(39) 


Like the KDM operator, it is self-adjoint, as we shall confirm. The extra term, involving 
reflection, arises because of the time integral in Eq. fl37)l . which is precisely how reflection 
effects arise in scattering theory in the usual S'-matrix expansion 

We also note that our operator bears a close comparison with the dwell time operator. 
This is the operator describing the time spent by a particle in a spatial region [0, L\, defined 
by 

poo 

T D = dt e* Ht P L e~^ H \ (40) 


where H = p 2 /2m and Pl = f 0 dx\x)(x\ is a projector on the region [0, L\ 32], The dwell 
time operator commutes with the Hamiltonian so is unaffected by the Pauli theorem and 
indeed it is self-adjoint. We would expect self-adjointness on physical grounds. In a mea¬ 
surement model for dwell time (see for example Ref. 24]), one might consider a clock model 
in which the Hamiltonian (which we take to be self-adjoint) includes a coupling between the 
clock and the operator Pi and the clock variables would then “see” the expression Eq. (140H 
in an S'-matrix expression, analogous to Eq. ([32]h 

The dwell time operator may also be written more explicitly as 

mL 


T d = 


\P\ 


<- L sin {pm R 
( PL/h ) K 


(41) 
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in which one can clearly see reflection effects. Classically, the dwell time can be written as 
the difference between two arrival times, at x = L and x = 0 


mL 

\p\ 


m(x — L ) 

\p\ 


mx 



(42) 


but not in general in the quantum case Eq. fldTjh However, it is true both for large momenta, 
pi h, where the quantum case essentially coincides with the classical one as one might 
expect, but more interestingly also for low momenta, pL h, where one can see that our 
new arrival time operator satisfies 


T, 


D 


mL 

\P\ 


1 + R 


= e* Lp fe~* Lp -T. 


(43) 


Hence the reflection term occurring in our new arrival time operator is closely related to a 
similar term appearing in the more familiar dwell time operator. 

Returning to our new time operator, we consider its covariance properties. Note that it 
has the same commutation relation with the Hamiltonian as the KDM operator, namely 

[H, T] = ih e(p) (44) 

since the extra term in Eq. (j39j) commutes with H. Applying the unitary Hamiltonian 
evolution operator to the eigenvalue equation, 

T\4>r) = r\(/>r) (45) 


we find that 

T (e~^ Ht \ 0 T )) = (r + e(p)t) (e~^ Ht \ 0 T )) . (46) 

For the KDM case, we can solve this equation for e~^ Ht \<fi T ) by splitting into positive and 
negative momentum states, thereby finding that the positive momentum states are shifted 
forwards in time and the negative momentum states backwards (as one can see also in the 
explicit eigenvalues, Eq. (|6|) ). so there is covariance in the positive and momentum sectors 
separately. However, for our new operator T, the presence of the reflection term mixes up 
the positive and negative momentum sectors and it can be shown that the spectrum does not 
transform in any simple way under unitary time shifts. This difference between T and the 
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KDM operator may also be seen by considering their commutation with the pseudo-energy 
£ = p\p\/2m. We have 

[if KDM ]=ih (47) 

which implies the simple shift properties stated above. However, the commutator with T 
has an extra term, 

K, T] = ift (l + ijf) (48) 

from which no simple covariance properties follow, as one can show. These results mean 
that we do not expect any obvious covariance properties in the spectrum of T. However, 
since the extra term is clearly quantum-mechanical in nature, we would expect covariance 
behaviour similar to the KDM case in the semiclassical regime. 


VII. SPECTRUM AND ARRIVAL TIME PROBABILITIES 


Consider now the spectrum of the arrival time operator Eq.(j38j). The eigenvalue equation 
in momentum space is 


rnh 


d_l_ 

dp \p\ 


1 + R ) + 


IpI 


1 + R 


d_ 

dp 


MP) = T<t>r(v)- 


(49) 


The eigenstates can be written as the sum of their symmetric and anti-symmetric parts 
4> t (p) — 4>t(p) + 0r (p) and the eigenvalue equation reduces to a coupled system of first order 
equations 


d_ 

dp 



mh 


(50) 


and 




(51) 


Solving for the antisymmetric part of the state we obtain the second order differential equa¬ 
tion 


^r(p) - + (^)V^(p) = 0. (52) 

It may be shown that this equation has two linearly independent solutions in terms of Bessel 
functions, an antisymmetric one, p|p|u/3 and a symmetric one |p|§ J_3 which 

is irrelevant so is dropped. Inserting this antisymmetric solution into Eq. (l5T]) and using 
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properties of the Bessel functions 33] we obtain (f>f(p) and the full normalized solution is 
then found to be 


1 

T 2 


Mp) = l |p|lj -> © +tpWij i ©)' 


(53) 


(Here, a possible overall factor of i is hxed by noting that the solution must satisfy 
4>*{p) = These eigenstates may be shown, at some length, to be orthonormal and 

complete and thus T is a self-adjoint operator. Using the asymptotic forms of the Bessel 
functions, we obtain approximations for the eigenstates in the large and low momentum 
regimes respectively, 


-ip) 


e-*‘ w (*h exp ( <f(p) £) £»1 


2rfU - 


\P I 


p 2 r 

2mh 


(54) 


« 1. 


k 4 ) (mh) ^ 

As anticipated, the eigenstates do not shift in a simple way under time translations, except 
in the large momentum regime, where they have the same covariance properties as the KDM 
states. 

The arrival time probability following from our new operator for a time interval [ri, t 2 ] is 
now given by 

r T 2 


J>(ti,t 2 )= / (ir |{V#r)| . 


(55) 

J T1 

We compare this directly with the measurement model, Eq. (l24l) . so we set T\ = 0 and consider 
the interval [ 0 , 72 ]. If t 2 is large compared to the energy time mh/p 2 , the time integral in 
Eq. fl55l) takes contributions from both large and small momentum regimes in Eq. (15411 . It is 
easy to verify that the large momentum regime will give the dominant contribution. The 
eigenstates in that regime are of the form (|6j) and thus the probability distribution is similar 
to that for Trdm, he. a variant of the Kijowski distribution 34 J]. Note, however, there 
is a difference in that there is a phase depending on e(p) which implies that the terms 
in the distribution representing interference between positive and negative momenta will 
have a different phase to that in the distribution derived from T^dm- This difference is 
not surprising given the presence of a reflection term in our operator. We thus get broad 
agreement with Eq. (j24ll . 

If the time t 2 is short compared to the energy time in Eq. (|55l) . the low momentum regime 
applies. The eigenstates are of the form (p T (p ) ~ C\p\, where C is read off from Eq. (15411 and 


16 









we easily obtain 


(56) 


|{V#» 


7r 


r 2 


2 (r (|)) 2 w-2^2 


(V#|<K£)|p||^), 


which is the kinetic energy density, Eqs. (j22i) . (1291) . the physically expected result in this 
regime. The numerical prefactors are not exactly the same as in Eq. (l29l) - they differ by 
about 20 percent but the overall factors in Eq. (156]) are fixed by the completeness relation of 
the eigenstates 0 T (p). There is no reason to expect perfect agreement since the two formulae 
have different origins. So again we get broad agreement with the measurement model. 

ffence our new arrival time operator Eq. (1381) gives the physically expected results for 
both large and small momentum regimes, and in particular, gives the behaviour anticipated 
by measurement models in the low momentum regime, behaviour not captured by any of 
the standard arrival time operators described in Section II. 

However, it is important to note that approximate agreement of the probability Eq. (l55|) 
with the low momentum arrival time formula Eq. (I22[) is achieved only for an interval [ 0 , 72 ] 
with T 2 small. This is all that is required since we are comparing with the formula derived 
from measurement, Eq. (j24l) . which concerns a time interval close to zero. The initial time 
r = 0 has a distinguished status in this model since the state is prepared then (and in 
the stopwatch model it is the initial time at which stopwatch and particle are in a factored 
state). Hence this feature is physically expected. It is also mathematically consistent since, 
as noted, the spectrum of our operator does not have any obvious covariance properties 
under time translation so there is no inconsistency in the model possessing a distinguished 
point in time. For more general time intervals [ 71 , 72 ] we get the low momentum result 
only for time intervals where both end-points are close to r = 0. For times Ti,T 2 which are 
very close, but large, the probability Eq. (l55|) is still an integral over a range of time of a 
distribution similar to the Kijowski distribution, not the low momentum formula. 

These features mean that the regimes in this model are not in general defined by the 
size of the interval |i 2 — Ti | compared to the energy time, as one might have thought from 
some measurement models, except for the special case of the interval [0,r 2 ]. To obtain a 
result of this form would require a more elaborate operator containing an explicit parameter 
modeling the precision of the measurement. This is of interest to construct and will be 
pursued elsewhere. 

Furthermore, as indicated earlier, there are two different motivations for looking at the 
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new time operator and the second motivation is the different question of exploring the con¬ 
sequences of taking an alternative but equivalent classical starting point, without necessarily 
making any claims about a connection to measurement models. In that case, our new opera¬ 
tor shows that a different classical starting point leads to a probability distribution |(^|0 T )| 2 
which is a variant of the Kijowski distribution for large momenta. At small momenta, the 
alternative starting point yields improved behaviour in the eigenstates, since the states go 
like \p\ instead of . This means that arrival time states obtained by superposing small 
ranges of r have finite (Aa;) 2 , and hence better localization properties than the four standard 
results discussed in Section II. So interpreted purely as an alternative classical starting point 
our new operator gives useful and physically reasonable results. 


VIII. SUMMARY AND CONCLUSION 

In the construction of arrival time operators, different classical starting points lead to 
inequivalent quantum operators with different physical predictions. The Aharonov-Bohm 
operator and its variants do not capture the expected physical behaviour in the low momen¬ 
tum regime, suggesting that a new starting point is called for. In this paper we constructed 
an arrival time operator taking a different starting point based on the current, inspired 
by measurement models and classically equivalent to — mx/\p\, thus obtaining an arrival 
time probability giving the expected physical behaviour in both large and small momentum 
regimes. In addition, we have shed some light on the intuitve origins of self-adjointness or 
its absence in both arrival times and the momentum operator on the half-line - it relates to 
the imprecision arising due to reflection when x is restricted to a half-infinite range. 
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